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1 Introduction 

There are numerous applications of integral geometry to various problem of image recon- 
struction, in particular to X-ray, emission, magnetic resonance, wave, acoustic, thermo 
and photoacoustic, scattering, Doppler tomographies as well as to radar technique and 
texture analysis. We propose here a general method of analytic inversion formulas for 
various transformations in integral geometry. Several examples for the case of integral 
transforms in a plane domain are given. 

2 Curves and integrals 

Let X and S be smooth 2-manifolds and $ be a smooth real function defined in A x E 
such that da;$ 7^ 0. For any point ex e E the set F (a) = {x G A; $ (x, a) = 0} is a 
smooth curve in A. We call $ generating function of the curve family {F (a) , a G S}. 

Let dS be the Riemannian area form; we define a Funk-Radon transform generated 
by $ by 

M^f {a) = [5 ($ (x; a)) fdS = I j^^ 

for continuous functions / compactly supported in A. The quotient fdS/d^^ denotes an 
arbitrary 1-form q such that d^^ A q = fdS. It is defined up to a term hdx^ where h is 
a continuous function. An orientation of a curve F (a) is defined by means of the form 
dx^ and the integral of q over F (cr) is uniquely defined. Let g be a Riemannian metric 
in a manifold A; the form dgS = ^^/gJdx) is the Riemannian line element. We have then 



M^f (a) 



F(a) |Vg$(x;(T) 
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where |Vga| = ^/g (da) is the Riemannian gradient of a function a. Suppose that the 
gradient factorizes through X and E that is |Vg$ (x, a) \ = m{x) fi (a) for some positive 
continuous functions m in X and /i in E. It follows that data of the Funk- Radon transform 
is equivalent to data of Riemannian curve integrals 

RJia)= [ /dgS,aGS 

JF{a) 

since Rf (a) = n {a) M$ (m/) (a) . The reconstruction problem of a function / from 
Riemannian integrals Rgf is then reduced to inversion of the operator M$/. 



3 Main theorem 

We assume now that S = M x S*^ and a generating function is linear in the first argument: 
^ {x; X, (p) = X + ip {x, (p) , X E M., ip E = {0 < Lp < 27r}. We call this function regular 
if (i) dxi' A dxi''^ 7^ in X X S*^ and (ii) there are no conjugated points, that is the 
equations i^, f)=4' {Vi and ip'^ (a;, Lp) = ip'^ (x, y?) are fulfilled for no x 7^ ?/ G X, 
ip G S\ 

Theorem 3.1 Let ^ {x; X,(p) = X + ip {x, (p) be a regular generating function in X x T, 
such that the function ip is analytic in X and the principal value integral 

N{x,y) = {P) r 2=Re / , (1) 

vanishes for any x ^ y & X. Then the reconstruction formula 

27rJo \Wgip{x,ip)\ 

holds for an arbitrary function f E L2 (-^)comp- integral ([^ converges in L2 (X)^^^. 
Remark 1. A more invariant form of ([2]) is a reconstruction of the densityt. 

Here the quotient dS/ D (x) depends only of the conformal class of the Riemannian metric 
g- 

Remark 2. G. Beylkin studied "the generalized Radon transform" in M" [3] which 
coincides with our Funk-Radon transform M$ in a Euclidean space (see also the last 
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Section). He constructed a Fourier integral operator which serves as parametrix for this 
transform. Beylkin's parametrix gives a high frequency approximation to a solution but 
not an exact inversion due to limitation of the method of FIO. The explicit form of the 
nucleus N {x,y) does not appear in Beylkin's approach but is important in our method. 
Note that if the nucleus does not vanish it is still a smooth function in the complement 
to the diagonal, since condition (ii) keeps the number of real zeros of {x, iVi 
constant. Then the right-hand side of ([2]) is a parametrix. 

Proof of Theorem. To simplify our arguments we assume that X is an open set in 
and g is the Euclidean metric, the area form is denoted dx. Let g be a 1-form in X 
such that dip A q = dx; we have Mf (A, (p) = J^^_^^^Q /<?• For an arbitrary x G X and any 
function / that vanishes in a neighborhood of x we find 

M./,A,.)dAd._ 



s 



$2(x; A,(^) 

27r 1 \ 

/ (y) g {y) A d^ (y, ^) 



X 



iipix,(fi) -ip{y,Lp)f 
Nix,y) f {y)dy 

X 

since dA = —dip as X + ip {y,ip) = and dA A g = —dy by definition of q. It follows that 
the function N is the off-diagonal nucleus of / and it vanishes since of the assumption. 
Therefore the nucleus / (x, y) of the operator / (see (j5]) is supported in the diagonal and 
by Lemma 13.31 below we have / (x, y) = a (x) (y) for some continuous function a. To 
calculate this function at a point xq G X we apply / to a density he that is equal to dx 
in e- neighborhood X^ of xq and /i^ = in the complement: 

Ihs (xo) = [ dx [ —— — — r-2 = /" ^ 



'X, Jo {ip{x,(p) -ip{xo,ip)) Jo ix, (V^(x,(/?) - V^(xo,</?)) 

Lemma 3.2 We have for arbitrary (p & S^, Xq € X and small e 

1^ / \ ■ f dx 2tt ^ , , 

D (xo, <p> = 2 = 2 + ^ 

JxAi^i^^v) -i^i^o,v)) \Vip{xo,ip)\ 

where O (e) < Ce where C does not depend of ip. 

Integrating over (y9 yields the equation D (xq) = — J {Vip (xq, 'p)\~'^ dip + O (e) . This 
implies since D (xq) does not depend of e. 

Proof of Lemma. Replacing ip by ip/ \ Vip (xq, v^)] we can assume that \Vip (xq, (p) \ = 
1. Shift the origin to the point Xq and suppose first that ip {x,ip) is linear in x. By an 
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orthogonal transformation we can assume that ip {x, (p) = xi. In this case we have 



D {xq, If) = Re 
- -2 



dx 



d,T2 



= —Re / dx 



£ xi + iO 
27r 



■\J e^—x. 



2_~2 



In the general case we can write -0^ (x, (p) = x\— p [x) , where p{x) = O {^xf^ is analytic 
in 3£-neighborhood of the origin if e is sufficiently small. We have 



11 . V- P' 



™2n+2 

n=l 



and 



dx 



Ix, 



dx 



+ a (xi + iO, X2) dx 



(4) 



Ix, (V' + «o)^ Jx. {xi + ioy 

The first integral in the right-hand side is equal to 27r and shall show the second integral 
is small. The field 

X2 



v{x) ^ \x2\ + y e"^ - \x\ , 



F2 



-Xi 



is defined and continuous for x e X^, 0:2 7^ 0. For any t, < i < 1 we define 2-chains 

f . . X + itv (x) , , 

Vt,±^ iz^ z {x, t) = — , ±a;2 > 0, \x\ = £ 



W, 



Xi + itd — xf, istxi 1 , < £, — 1 < s < 1 



oriented by the orientation of X^. The boundary of the chain Vt — + Vt- + Wt is the 
image of dX^ under the map x ^ z {x,t) . The function a has a holomorphic continuation 
to the domain Z+ = {2; G C^; \z\ < 3e, yi > 0}. The form adz is holomorphic in Z+ hence 

1 /2 

dcr A d^ = 0. Consider an open 3-chain Zs — ^t=s^t C Z where 0<5<l/2. Itis 
contained in Z^ since yi ^ Im zi — t ^X2\ + \J ^'^ — > and \v\ < 2e. By Stokes' 



we have 



adz 



dcr A d2; = 



The boundary of is equal to V1/2 — + 5* + T where 

X + itv {x) 



S^lz^ 



T^lz^ 



{e, iste) 



{stf 



, \x\^e, 5 <t< 1/2 j , 
,e^±l,5 <t<l/2,-l<s<l 
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Therefore 



adz 



Vi/2 



Vs 



adz = 



S+T 



For any point z E S we have 

zl + Z2 



) = (1-^^) 



\X\ 



j-2 I |2 



) 



since — = and = on dX^. The same equation holds in T. It follows that 

the chain S' + T is contained in the complex algebraic curve z'^ = 6^. This implies that the 
integral of the holomorphic form adz over S + T vanishes. This yields for an arbitrary 
5 > 

adz = adz = J2 

Vs JVy2 n=l -^^1/2 ^1 

The left-hand side tends to the second term of (jlj) as 5 — )■ 0. The right-hand side can be 
estimated from above by the sum 



E 

n=l 



p^dz 



Vl/2 ^1 



2n+2 



< 



I d5VM"max|zi| 



-2n-2 



where 



M. = maxlpl < Ce^ 



and dS is the Euclidean area element. The area of the chain Vi is estimated by const 
since |Vf | is bounded in X^. For we have any point z E Vi 



This yields 



adz 



1/2 



< 



(3M, 



Ce 



e 



2{n+l) 



if £ < C ^. This completes proofs of Lemma [3.21 and Theorem 13. 1[ 
Lemma 3.3 The integral transform 



If (x) = (P) 



27r 



M^fiX,^) 



is a continuous operator L2 (-^)comp ^ ^2 (-^)ioc • 



dXdip 



(5) 
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Proof of Lemma. 

where 



The condition (i) can be written in the form J^-,^ ($) 7^ in F 



Jx-xv i^) = det 





92$ 






dxidX 


dxidip 


dxi 








9$ 




dx2dX 


dx2dip 


8x2 


v 


d<l> 


a* 




dX 


dip 






(6) 



xi,X2 are coordinates in X. According to [4j, Theorem 25.3.1 the map : L2 (-^)comp ~^ 
L2 (S)jqj, is a Fourier integral operator of order —1/2 since the corresponding Lagrange 
manifold is locally a graph of a canonical transformation (see details in [TU] ; the bundles 
^1/2 (^x) and fi^/^ (S) are trivial in our case). The image of this operator is contained in 
L2 (S)^Qjj^p since for an arbitrary function / G -L2 (-^)comp support of M$/ is contained 
in 9 = TTsvr^^^ (supp /) where ttx and tts are projections of F to X and S respectively. 
The projection ttx is obviously proper. The integral transform 

h (A, (f) dXdip 



Sh 



x] 



$2 {x;X,ip) 

can be written by means of a Fourier integral operator 



Sh< 



x] 



- / exp {it^ {x, X,ip)) \t\ h {X,(p) dtdXdip 
^ Jt, Jr 



with the phase function $. The corresponding Lagrange manifold coincides with that of 
M$ with source and target spaces interchanged. The order of S equals 1/2 due to the 
factor \t\. It follows that the composition I = SM^p is well defined as a PDO of order 
0. ► 



4 Integrals of rational trigonometric functions 

We study now the condition N (x, y) = 0. A function 

k 

t {(p) = ttm COS rmp + bm sin rmp 

m=0 

is called trigonometric polynomial of order k if la^f + \bkf > 0. Any trigonometric 
polynomial is 27r-periodic and is well-defined and holomorphic in the cylinder C/27rZ. It 
always has 2k zeros in the cylinder. If a polynomial is real the number of real zeros is 
even. 

Lemma 4.1 Let t (9?) , s (if) be real trigonometric polynomials such that deg s < n deg t 
for a natural n and all zeros of t are real and simple. Then 
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Proof. The function r {() = s {C)t~"' {() is meromorphic for ( = ip + ir E C/27rZ 
and has no nonreal poles since of the assumption. Fix a small positive e we choose a real 
continuous function A = r {(f) defined on the circle < ip < 27i that vanishes except for 
e-neighborhood of the zero set of t and 

T (if) = sgn t' (a) sje'^-iyif- af 

in e-neighborhood of each zero a of t. We have t{a + it (a) 0) = t {a) + zO for any zero 
a and 

Jo {t{^)+tor Jo t-{c) '^^^ ^ 

This yields 

f _i(^d^+ r_i(^d^= / 4id^ 

i„ (i(v) + jO)" ^ i„ (((v')-«0)" irur«"(C) 

where F = = ^ (ip)} . We can write F U f = F+ U F_ where T± = {( = ip ± i\T {(p)\} . 
Finally 

since t (v^ + «t) 7^ for r 7^ and the function tends to zero at infinity. ► 
The following formula can be useful for calculation of D (x) : 

Lemma 4.2 If t and s are real trigonometric polynomials, deg s < degt = k and t has 
no real zeros, then 

where the sum is taken over k zeros of t with positive ( or negative ) imaginary part. 

For a proof we apply the Residue theorem for the form s {ip) dip/t [ip) in the upper 
half- cylinder C+/27rZ. 



5 Radon's and Funk's reconstructions 

We apply Theorem 13.11 to recover some known and unknown formulas for regular curve 
families. 

Radon's formula. Take a generating function $ (x; X,uj) = X — {x, e (</))) in ]R2 X S 
where e (ip) = (cosy), simp). The classical formula of Radon- John 



coincides in this case with ([2]) since integral data g = R f are equal to Mf (A, ip) because 
of iVipl = 1. The equation N{x,y) = immediately follows from Lemma l4.ll The 
coefficient is reduced to — l/27r^ due to symmetry g (—A, ip + n) = g {\,ip). 

Punk's formula provides reconstruction of an even function / in the unit sphere S'^ 
from integrals of / over the family S of big circles. Take the unit hemisphere X = {x & 
E^, \x\ = 1, xo > 0} and a generating function $ (x; A, = A + ■?/', ip = {Vi ^ (v^)) where 
Hi = xi/xo,y2 = X2/X0 defined in X x S. By Lemma [4.11 we have N{x,x') = for any 
X 7^ x' G X. Let g be the standard spherical metric in X. We have 



Xq 



dif 2 27rxo 



Xn / Z J . = , = 27TXr^ 



iVg^i' Vo i + y^cosv yrT7 

By Theorem 13.11 for any function / € L2 (X) 

iTT^xlJ J (A+(y,e(v)»' 

Choose coordinates 6, — 7r/2 < 6 < 7r/2, in the dual sphere S so that a = (sin 6, cos^^ cos cos 6* sin ip) G 
E. We have A = tan 6', X + {y,e = {a, x) /xq cosO and |VgV'|~"'^ = Xq cos6'. 



M/(A,^) 



/ tItTi = cos 6* /" xo/dgS = cos9 R (xo/) (o") 



We apply (jH]) to the function h (x) = Xq/ (x) and rearrange this formula as follows 
1_ i?g/((T)cosgdgdy? - J_ f ^sf (^) 



where da = cos^d^d<y9 is the area form in the hemisphere 5"+. This formula coincides 
with original Funk's result P after a partial integration. 



6 Geodesic transform in Lobachevski plane 

Take a generating function ^ {x; X, (p) = X + ^jj, = —2 (|x|^ + l) ^ (^5 ^ (v^)) 5 ^1 < 
A < 1 in the unit disc X = D. The curves F (A, y?) are geodesies in the Poincare model 

1 — |x| j ds is the hyperbolic metric. We have 
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and 



Mf (A, V.) = (1 - Xr'^' [ / (x) d. = (1 - Xy'^' [ ^-^f (x) d,., 



where ds, dgS means the Euchdean and hyperbohc hne element respectively. Further 

2nJo \Vij{x,^)f Stt^ ''Wo I - cos^ ^ 4 (l - 

since the integral in the right-hand side is equal to 2n (l — \yf) ^ , y = 2 \x\ (^1 — \xf) 
(e.g. Lemma [4. 2p . By Theorem 13 . 1 1 this yields 

We apply this equation to the function /* (x) = 4 (l — |x|^) / (x) and get 

. . ^ _^ _ , f Rf{X,^){l-Xr'^'dXd^ 

4vr2 ^ ' ' Wo J-i ((l + |x|^)A-2(x,e(v,)))' 

where R f {X,ip) = J fd^s and the integral is taken over a geodesic F (A, ip) in Lobachevski 
plane. This formula is equivalent to Helgason's theorem [6]. 
Applying a map 

y' 

X = — — , y' = (z/i,z/2) 
1 + Z/o 

we transform the Poincare model to the Lorentz model in the " upper" sheet hyperboloid 

Q+ = {{yo, y) e E\ yl = l + yl + yl > 0} 

with the metric induced from the Euclidean metric in . For any central plane P in 
the hyperbola 7 = P fl is a geodesic curve in Q^. Vice versa, the image in (5+ 
of an arbitrary geodesic circle F (A, = {A = 2 (l + |x|^) (x,e (</)))} is contained 
in the plane P = {y; Xyo — (e {^p) ,y') = 0}. This plane is orthogonal to a vector ^ = 
(1 — A^) ^^'^ {X, —e {(p)) which is contained in the one-sheet hyperboloid = G 
E^,C,Q — ^1 — ^2 = Therefore we use notation 7(^) = F{X,<f) = P H for a 

geodesic in and the hyperboloid parametrizes all the geodesies in the Lorentz 
model. We can write $ (x; X,(p) = X — yQ^ {y, e ((/?)) = (1 — A^)"^^^ y) and have 
1 — |x|^ = 2 (1 + 2/0)^^ 5 1 + = 2?/o (1 + yoY^ ■ The equation (|TT]) is now read as 
follows 

1 f f Rf (0) {I - xT'^' dxd^ 



Expressing the projective area form 

u (0 = eod^id^ + 6de2deo + ^d^odei 
in terms of coordinates A and ip we get u (0 = (1 - A^)"^/^ dAdy?. Finally 

j^y^_ 1 /■ i?/ (7 (0)^(0 



which makes complete similarity with ( ITOl) . 

Remark. A formula of this form appears in [7], p. 96 however with misprinted nu- 
merical coefficient. 



7 Equidistant curves in Lobachevski plane 

Let X = D he the open unit disc and F = {F{X,ip)} be the family of arcs in D 
connecting two opposite points e{(f — 7r/2) and e{(f + 7r/2) on dD where A = tana;, uj is 
the angular measure of F (A, ip) . For any ip the arcs F (A, ip) and F (A', ip) are equidistant 
with distance 2 |A — A'| . This family has generating function 

$(x;A,y.) = A + ^, ^{x,^) = -^-p^ML 

1 — \x\ 

Check that this function fulfils the conditions of Theorem 13. 1[ A proof of regular- 
ity is a routine. Further we have ■ilj{x,ip) — ip{y,ip) = {y' — x' , e (ip)) where x' = 

1 — X, y' = [1 — \y\ ) y. This is a first order trigonometric polynomial with 

two real roots which implies N {x,y) = as x ^ y. We have 

(x, ip)f = i— 2 (1 + cos^ ip) 

{l-\x\ ) 

where z = 2 \x\ (^1 — ^ and 



Stt Jq l + z^cosV 4 

Corollary 7.1 For any function f with compact support in the unit disc a reconstruction 
from M^f is given by 

fix) ^ -ii±^ r r ^Mi^ ^dAd, 
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Because of |V$(x;A,v9)| = 2 (l — ^ cos u we can express the Funk-Radon 
operator M in terms of Euclidean arc integrals 

Mf (A, if) = ^Rf* (A, ^) , A = tana;, /* (x) = (l - \x\') f (x) 

This yields 

l + \xf i?/(a;,(^)cosa;da;d(p 



fix) 



27r2 



J-tt/2 ((l — |a;| ) sinw — 2 {x, e cosw) 



This formula was stated in ^HJ by a different method. 



8 Photoacoustic geometries 

Elliptical source curve. D. Finch and Rakesh [S] gave a formula of type for 
reconstruction for the family of spheres centered on a sphere. Another reconstruction 
formula was proposed by L. Kunyansky [9]; it looks different however can be reduced 
to a form close to (|2]) after one-fold integration. We show that a reconstruction can be 
done for ellipse and more general source curves by means of Theorem 13.11 F. Natterer 
applied quite different method for reconstruction from the family of spheres centered in 
an ellipsoid [12]. 

Take first a generating function ^ = X + ip where 

ip {x, (fi) = \x — e {(f)f , e {(f) = (cosyjci, sin (^62) , 

where 61,62 are arbitrary positive numbers. The source curve s = e{ip) , ip G is an 
ellipse with half-axes ei and 62- 

Proposition 8.1 For any function f with support in the ellipse E = {x; (xi/ci)^ + 
(2^2/62)^ < 1} a reconstruction is given by 

f (^) = -2 — / / 7 , TT^^^^*^ 



{x-\x-e{ip)\ y 



J J (r'^ — \x — e {(p)f)'^ 
where e = {ci, 62) and A = since Mf (A, ip) = R f (r, ip) /2r, R f (r, 9?) = /p(^_^) fds. 
Proof. We have 

^ (x, ^)-ip (y, ip) = \x~e {ip)f -\y-e{ip)f = 2{y-x,e (if)) + - \yf 



2 \\y — cos {(f — 6) + \xf — \yf (12) 
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where 9 = arg {y — x) and \\z\\^ = ((2:161)^ + (-2262)^)^''^ • This norm and the dual norm 
Iklle ~ ((-^i/ci)^ + (-22/62)^)^^^ fulfil the inequality 

\\xf - \yf\ < \yi -xil \yi + xi\ + |?/2 -3^2! \y2 + X2\ 
< \\y - x\\e\\y + ^Wl < 2||y-a;||g 

since \\y + < 2 for x, y belonging to the support of /. This implies that the trigono- 
metric polynomial ( JT2l) has only simple real zeros. By Lemma 14.11 [x, y) = for 
arbitrary x y E E. We have = 2\x — e {{p)\ and 



D(x) 



dip 



7f J + |e|^ — 2 |x| |e| cos (v9 — 6*) 4 



\x\ 



and our statement follows. ► 

A similar reconstruction from spherical means in can be done by means of Theorem 

[im 



9 Isofocal hyperbolas and parabolas 

Hyperbolas. The equation A = |x| — exi,6 > 1 defines a fold of the hyperbola 



axi H — ) — Xn = a = Ve^ — 1 
a J a'' 

with a focus at the origin. The function $ (x; A, ip) = X+ip (x, (p) , ip {x,ip) = e {x, e (yj)) — 
|x| generates the family of all one-fold hyperbolae with focuses in the origin. The first 
order trigonometric polynomial 

^ (x, (p)-^ (y, (p) = e {x -y,e {(p)) - \x\ + \y\ 

has two real zeros if x y. since ||x| — \y\\ < e \x ~ y \ . By Lemma HTT] the nucleus (x, y) 
vanishes for y ^ x. We have {Vipf = 1 + e"^ — 2e \x\~^ (x, e) and 

D (x) 



2n Jo l + 52_2e|xr^(x,e(^)) £^ - I 
(Lemma 112]) • 

Corollary 9.1 For any smooth function f with compact support in the equation holds 

f^^^^ e'-l r Mf (A, ip) dXdip 



4vr2 Jo i_oo (|x| -£(x,e((^)))^ 
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Parabolas. A parabola with a focus at the origin can be given by 2pxi +x\—p^ = Q 
where p is a positive parameter. Set A = p^/"^ and write this equation in the form 
(xi — A^)^ = \x\^ which is equivalent to \x\ + Xi = A^ since Xi < A^. Take a function 
$ (x; A, v?) = A + V (a;, ip) in X = R2\{0} where 

if) (x, if) = — a/|x| + (x, e (ifi)) = — (2 Ixl)"*"^^ cos {(p — 6) /2, 9 = arg x 

This function generates all parabolas with focus at the origin. For any x ^ y & M? 
the trigonometric polynomial (x, (p) — ip {y, (p) is of order 1 and zero constant term. 
Therefore it has two real roots and by Lemma WA\ the function (x, y) vanishes iiy ^ x. 
We have 

= Mf (A, cp)= [ (2 \x\f' f (x) ds, D{x) = 2 \x\ 

Corollary 9.2 For any L2-function f with compact support in M^\{0} a reconstruction 
is given by 

^""^ 1^1-^0 (^X-^\x\ + {x,e{^))y 

^ 1 r Rfix,ip)dxdip 

4^^H'2\x\f^ Jo Jo (^X-^\x\ + (x,e((^)))' 
In [S] Fourier coefficients of / are reconstructed from Mf by Cormack's method. 



10 Cormack's curves 

For a > an a-curve is given in a plane by a equation 

A = r" cos {a {9 - , \9 - ^\ < 7r/2a 

where x = re {9) and A, ip are parameters Suppose that a = A; is natural and take the 
family of curves generated by a function $ = A + (x, , (x, = —r^ cos {a9 — ip>) . 
Note that in the case k = 2 the curve ijj (x, (p) = —A is a hyperbola. The function ijj 
always admits a symmetry group acting by rotation of the plane by the angle 27r//c 
about the origin. Consider a quotient manifold X = M^\{0}/Zfc; lifting of $ to X is a 

regular generating function. We have ip (x, ip) = Re ^exp {—i<p) (xi + 1x2)^^ and 
ilj{x,(p) -ip {y, ip) = Re exp - i<p ( (xi + ixs)*" - (yi + ^2)'' ) 
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is a first order trigonometric polynomial with two real zeros for arbitrary x ^ y. This 
implies (x, y) = 0. Further we have 

\V^{x,^)f = e\xf'-\ Dix) = —^^, Mf{X,^) = l [ 

\x\ J\=tp{x,<p) \x\ 

and the reconstruction formula reads: for any Z^-invariant function / with compact 
support in M^\{0} : 



fix) 



47r2 



' — Reexp (—i^) (xi + ix^f" 
The family of /3-curves A = cos f3 {6 — (f) with natural (3 is studied in the same way. 



11 Higher dimensions 

Let $ be a generating function defined in X x S where X C M" is a smooth manifold of 
dimension n with a Riemannian metric g, S = M x S""^ and ^ {x; X,uj) = X + ip {x, u), 
A G M, w G S*"^^. We say that a generating function $ is regular ifan+lxn + 1 
determinant like ([6]) does not vanish in F = {$ = 0} and conjugated points are absent, 
see [in] for details. The Funk-Radon transform is defined by 

fdV 



M^f{X,ip) 



F(A,0) 



where dV means the Riemannian volume form and F (A, = {x : $ {x; A, = 0} . The 
result of Sec. 3 is generalized as follows. Define a function 

N{x,y)=Re — ^ ^ n , -n^" '^ ^ ^ 

J sn-1 [i) (x, uj)-i) [y, uj) ± lO) 

where do; the Euclidean volume form in S"""^ and 

1 f du 



D x) 



Theorem 11.1 If a regular generating function $ is real analytic in X and N {x, y) = 
for any x ^ y E X then a reconstruction of type ([^ exists for an arbitrary function 
/ G L2 (X) with compact support. For even n it reads 



(27ri)"D(x) Jj^{^{x]X,uj) + iQf 
For odd n it is 

fix) = }_.^ ^ [ 6^--'H^ix;X,u))M^fiX,u)dXdu 



2{2mf-'D{x) Js 
The integrals converge in L2 (-^)ioc ■ 
A proof will be given elsewhere. 
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